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Abstract 

We report on results of a search for the decay mode Kl — > ir je + e~ conducted by 
the E162 experiment at KEK. We observed no events and set a 90% confidence level 
upper limit of Br{Ki — > 7r°7e + e _ ) < 7.1 x 10~ 7 for its branching ratio. This is the 
first published experimental result on this decay mode. 

PACS: 13.20.Eb, 14.40. Aq 



1 Introduction 

Chiral perturbation theory (%PT) is a very powerful tool to describe vari- 
ous K decays in which long distance contributions are expected to dominate. 
For example, the decay mode K + — > 7r + 77 has been observed recently, and 
compared with [1]. The measured branching ratio and ir + momentum 

spectrum are found to be consistent with the predictions, after fitting one 
free parameter contained in the theory. The neutral counterpart, Kl — > 7r°77, 
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is another decay mode, where detailed studies have been performed. In this 
case, the lowest order 0(p 4 ) calculation [2] does not reproduce the measured 
branching ratio [3-5]. Extending the calculation to the next-to-leading order 
0(p 6 ) [6], and adding a vector meson contribution [7], the prediction is now in 
good agreement with the branching ratio as well as distinct M 77 spectrum. We 
note that an effective coupling constant (ay), a free parameter in the vector 
meson contribution, has been determined by the measurements with experi- 
mental errors [3,5]. The as yet unobserved decay mode Kl — > 7r°7e + e~[][] can 
provide another testing ground for xPT. Since theoretical ingredients are same 
for both Kl — > 7r°77 and Kl —> 7r°7e + e~ modes, a straightforward extension 
from Kl — > 7r°77 provides definite prediction for a branching ratio and M iee 
spectrum [8]. In particular, the branching ratio is calculated to be 2.3 x 10 -8 . 
Thus experimental study of this mode is important to test the theoretical 
framework of ^PT. 

Other interests in this decay mode stem from its close relationship to the 
mode Kl — > 7r°e + e~, which has been of much attention as a possible channel to 
observe direct CP-violation. First, Kl — > 7r°7e + e~ is expected to have a much 
larger branching ratio than Kl — > 7i°e + e~, and hence can be an experimental 
background in a soft photon region. Second, there also exists a CP-conserving 
amplitude in Kl — > 7r°e + e~ via two-photon intermediate states; this can be 
in principle determined from a detail analysis of Kl — > 7r°77 [9]. Further 
understanding of the Kl —>■ 7r°77 amplitude, which can be checked by the 
Kl — > 7r°7e + e~ mode, is thus essential. 

In this article, we report on an experimental search for the decay mode Kl — > 
7r°7e + e~ conducted with a proton synchrotron at High Energy Accelerator 
Research Organization (KEK). 



2 Experimental setup 

The data for the Kl — > 7c 0/ ~fe + e~ mode were recorded simultaneously with the 
experiment which has established a new decay mode Kl — > 7r + 7r~e + e - [10]. 
Since the experimental set-up was described already in Ref. [10,11], it is briefly 
mentioned here for convenience. 

The K L beam was produced by focusing 12-GeV/c primary protons onto 
a 60-mm-long copper target. Its divergence was ±4 mrad horizontally and 
±20 mrad vertically after a series of collimators embedded in magnets. The 
set-up started with a 4-m-long decay volume. It was followed by a charged 

4 An observation of 18 events has been claimed by the KTEV experiment in ICHEP 
98 workshop 
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particle spectrometer consisting of four sets of drift chambers and an ana- 
lyzing magnet with an average horizontal momentum kick of 136 MeV/c. A 
threshold Cherenkov counter (GC) with pure N 2 gas at 1 atm was placed in- 
side of the magnet gap to identify electrons. For the present decay mode, we 
obtained an average electron efficiency of 94% with a pion-rejection factor of 
350 by adjusting software cuts in the off-line analysis. A pure Csl electromag- 
netic calorimeter played a crucial role in this analysis. It was located at the 
downstream end of the spectrometer, and consisted of 540 crystal blocks, each 
being 70 mm by 70 mm in cross section and 300 mm (~ 16X ) in length. It was 
configured into two banks of 15 (horizontal) x 18 (vertical) matrix. Its energy 
and position resolutions were found to be approximately 3% and 7 mm for 
1-GeV electrons, respectively. The trigger for the present mode was designed 
to select events with > 2 electron-like tracks and > 3 cluster candidates in 
the calorimeter. It was formed with information from GC and Csl together 
with trigger scintillator hodoscopes interspersed between the chambers and 
calorimeter. 



3 General event selection 



In reconstructing Kl — > ir 'ye + e~ events, major backgrounds are expected 
come from both Kl — > 7r 7r^, and Kl — > 7r°7r°7r^ modes, where denotes 
the Dalitz decay n° — > e + e~7. Among them, Kl — > is an intrinsic back- 

ground and can not be removed. It is instead used as a normalization mode 
by positively identifying the Dalitz decay. The Kl — > n Tr it D mode may be- 
come background when two (or more) photons fuse into one in the calorimeter, 
and/or its vertex is reconstructed incorrectly to give false invariant masses. 
Care was taken in this analysis to enhance position resolution of decay vertex 
and purity of photon clusters. This effort was found useful also to reject back- 
grounds originating from external conversions and hadron (mostly neutron) 
interactions. 

In the actual off-line analysis, we first selected events containing two tracks 
with a common vertex in the beam region inside the decay volume. We then 
requested events to have > 5 clusters in the calorimeter. Here the cluster was 
defined as 3 x 3 Csl blocks around the local maximum with the total energy 
deposit greater than 200 MeV (± 3.5 nsec timing window). Particle species 
were then determined. A charged track which could project onto a Csl clus- 
ter was called a matched track. An electron (or positron) was identified as a 
matched track with 0.9 < E/p < 1.1, where E was an energy measured by the 
calorimeter and p was a momentum determined by the spectrometer, respec- 
tively. Also GC hits in the corresponding cells were requested. Clusters which 
did not match with any charged tracks were treated as photon candidates. 
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Event topology was then checked. We requested events to contain exactly one 
e + e~-pair and three photon candidates, consistent with the Kl — > 7r°7e + e~ 
topology. No additional activities, such as an extra track, GC hit, or cluster 
with an energy above ~60 MeV, was allowed in the detectoif^]. The probability 
of over-veto was estimated using Kl — > 7r°7r^ and Kl — ► e + e~ / y reconstructed 
events, and was found to be about 12% (common to both Kl — > 7t 0/ -fe + e~ and 
Kl — > 7i°7i% modes). 

Two quality cuts were imposed at this stage. One was a cluster shape cut; 
it examined mainly a cluster's transverse energy profile, and checked whether 
or not it was consistent with a single photon. We used a large sample of 
Kl — > 7r°7r + 7r~ reconstructed events to characterize actual photon showers, 
and found the cut efficiency to be 95% per a single photon cluster. The other 
quality cut was imposed on an e + e~-track opening angle (9 ee ). Larger opening 
angle resulted generally in better vertex position resolution, which in turn led 
to better invariant mass resolution and background rejection. We employed a 
Monte Carlo simulation to study background rejection power, especially for 
the Kl — > 7r 7r°7r£, mode, and determined to demand 9 ee > 20 mrad. 

As a final step in the general event selection, we imposed two loose kine- 
matical cuts to reduce sample size. They were M eery7 > 400 MeV/c 2 and 
9 2 < 100 mrad 2 , where 9 represents the angle of the reconstructed eeyy-f 
momentum with respect to the line connecting the target and vertex. 



4 Kinematical reconstruction 

Having selected candidate events, we scrutinized each event from the view- 
point of kinematical variables. First of all, we calculated 77 invariant mass for 
three possible combinations. The solid- line histogram shown in Fig.l is the dis- 
tribution for all combinations while the shaded one for the combination closest 
to M n o. We retained all combinations which satisfied |M 77 — M n o\ < 2.5er 7r o0, 
where a^o being the observed mass resolution of 5.1 MeV/c 2 . 

We then calculated the angle (9*) between reconstructed 7r° and ecy momen- 
tum vectors in the Kl rest frameQ. Events originating from Kl — * 7r°7e + e™ or 

5 The exception was bremsstrahlung photons. When track segment upstream of 
the magnet could be projected onto a neutral cluster with an energy below 400 MeV, 
we simply ignored this activity. 

6 This condition was determined in such a way that the number of final Kl — ► 

background events, estimated by a Monte Carlo simulation, became less 
than one inside the final signal box (see below). 

7 Actually this frame was obtained with the Lorenz boost, along the line connecting 
the target and vertex, defined by the velocity of Kl with an observed ee'yyy energy. 



4 



> 

CD 



CD 
> 
LU 



225 
200 
175 
150 
125 
100 
75 
50 
25 




all combinations 
closest to M 




0.09 0.1 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 

(GeV/c 2 ) 



Fig. 1. The invariant mass distribution of 77. The solid-line histogram shows the 
distribution for all 77 combinations while the shaded one for the combination closest 
to M^o. The arrows indicate the cut position. 

K L — > ir°7r% should satisfy cos^* = — 1 in this frame. The dots with error bars 
in Fig. 2 show the cos 6* distribution; a clear peak of events at cos 6* = — 1 can 
be seen. The histogram in Fig. 2 is Monte Carlo data for the Kl — > 7r°7r°7r£, 
mode, in which the flux was normalized by the observed K L — > -k°-k° d events 
(see below). To select signals, we requested events to satisfy a collinearity cut, 
cos6* < —0.98, as shown by the arrow in Fig.2. 

We now identify Kl — > h°t{ q d events. If more than one 77 combination within 
an event satisfied the M n o cut, we selected the one for which the quantity, 



M ee7 - M n o 



M~~ - ikLo 



became minimum. Here cr_o is the observed 7r mass resolution for the Dalitz 

D 

decay mode (see below for the actual value). 

Fig. 3 (a) shows a scatter plot of M ee777 vs 9 2 after selecting the 77 combination 
with Xp- To select K L — > 7r°7r^ events further, 9 2 < 30 mrad 2 was demanded. 
Fig. 3(b) shows a scatter plot of M ee777 vs M ee7 after this cut. A clear cluster 
of K L — > ti q ti d events can be seen in the expected region of M ee7 = M n o and 
M ee777 = Mk l - Fig.3(c) is the projection onto the M ee7 axis, and Fig. 3(d) 
is onto the M ee777 axis with a ix° D mass cut (see below). From these pro- 
jections, we found -k° d and K L mass resolutions to be ~ 4.5 MeV/c 2 and 
<jk l — 16 MeV/c 2 , respectively. Our final signal box, shown by the rectangle in 
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Fig. 2. The distribution of cos8* . The dots with error bars show the experimen- 
tal data, and the histogram is Monte Carlo data for Kl — > 7r°7r°7T£ ) , in which 
the flux was normalized by the observed Kl — ► TT°Tr% events. The cut position 
(cosO* < -0.98) is indicated by the arrow. 

Fig.3(b), was defined by \M eei — M^o\ < 3a^o and |M ee777 — Mk l \ < 3&k l - Af- 
ter all the cuts, 49 events remained. We estimated the number of Kl — > 
backgrounds in the signal region by a Monte Carlo simulation, and found to 
be less than one. 



We are now in a position to look for the Kl — > tt ^ye + e~ mode. First we 
rejected Kl — > ti ti° d events; if an event satisfied both |M 77 — M n o\ < 5a n o and 
|M ee7 — M n o\ < 5(7^ for any 7 combinations, then the event was discarded. 
Note that we employed the looser kinematical cut of 5cr 7r o(cr 7r o ) ) to exclude 
possible Kl — > events. Then we rejected Kl — > it Tt K° D backgrounds. In 
this case, an event containing 71% whose transverse momentum was consistent 
with Kl — > 7t°7t°7T£, (pt < 139 MeV/c) was excluded. The remaining events 
are shown in a scatter plot of M ee777 vs 9 2 in Fig. 4(a). We are left with 
no events inside our signal box defined by 9 2 < 30 mrad 2 and |M eerry — 
Mx L I < 3o"^ L Q. The background events still remaining in the low mass region 
of M ee777 < 460 MeV/c 2 were found to originate mostly from the K L — > 
7T°7T°7r^ mode. The projection of events with 9 2 < 30 mrad 2 onto the M ee777 
axis is shown in Fig.4(b). The solid-line histogram is for the data, the shaded 
one for Monte Carlo events (sum of Kl — > tt°tt^ and Kl — > ^tt ^) and 
the dotted one for Kl — > 7r°7e + e~ Monte Carlo events. From the Monte 

8 The actual value of a* K was 11.6 MeV/c 2 . This was the value of the Monte Carlo 
mass resolution (8 MeV/c 2 ) for Kl — ► 7r°7e + e _ times the ratio of the observed 
(16 ± 2.6 MeV/c 2 ) to Monte Carlo (11 MeV/c 2 ) resolutions for K L -» vr ^. 
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Fig. 3. (a) The M ee7Ty vs # 2 scatter plot of the Kl — > tt^tt^ candidate events, (b) 
The scatter plot of M ee777 vs M ee7 for the events with 6 2 < 30 mrad 2 . The box 
indicates the signal region, (c) The M ee7 projection of the events in (b). (d) The 
M ee7TY projection of the events in (b) with the ir° D mass cut (see text). 



Carlo simulation, whose flux was normalized by the observed Kl 
signals, we expected 1.1 backgrounds (0.45 from Kl 



7T 7T 



D 

7r° 'n° D and 0.66 from 



Kl — > 7r 7r 7T£,) to remain in the signal box. 



5 Results 



The branching ratio is calculated by 



Br(K L -> 7r°7e + e-) = 5r(fT L -> tt ^) 

A(7r°7e+e") ' r7(vr°7e+e-) ' ^(tt ^) 
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Fig. 4. (a) The scatter plot of M eerry vs 6 2 for Kl — > 7r°7e + e~ candidate events. 
There are no events left inside the signal box. (b) The M ee777 projection of the 
events with 9 2 < 30 mrad 2 for the data (solid-line), Monte Carlo (shaded, sum 
of Kl — > vr ^ and Kl — > 7r°7r 7r£,) and Kl — > Tr°je + e~ Monte Carlo (dotted, 
arbitrary scale). The arrows show the signal region. 



where A, 77 and N denote acceptance, efficiency and observed number of events, 
respectively. The detector acceptances depend on the matrix elements: we 
employed a theoretical model given by Ref. [8] for Kl — > 7r°7e + e _ and Kroll- 
Wada spectrum for the Dalitz decay [12]. The actual acceptances (including 
the hardware trigger conditions) were determined by Monte Carlo simulations 
and were found to be 2.34 x 1CT 4 for K L — > 7T°7e + e~ and 1.38 x 1CT 4 for 
K L — > 7r°7i£), respectively. The detector efficiency was found to be practically 
same for both modes. However, some of the quality and kinematical cuts 
caused efficiency differences. 

Since the opening angle of e + e~-tracks was different for the two modes, it 
caused efficiency difference in both the vertex reconstruction and the track 
opening angle cut. These efficiencies were found to be 70% (vertex) and 89% 
(9 ee ) for Kl — > 7r°7e + e~ mode while 62% and 81% for Kl — > ti ^ mode, 
respectively. Next, the collinearity (8*) cut produced 25% inefficiency for the 
signal mode while negligibly small loss for the normalization mode. The n° D 
mass cut was unique to K L — > tt ^; its efficiency was found to be 96% (3a 
cut). For the signal mode we rejected K L — * vr ^ events; it caused 9% ineffi- 
ciency at the specific ee'-f invariant mass region. We also rejected Kl — > 7r 7r 7T£, 
events with the tt% inclusive veto; it resulted in 7% inefficiency Combining 
other efficiencies together, we found the final acceptance and efficiency ratio 
A^Tr^/A^TT^e+e') to be 0.670. 



Inserting the known branching ratios [13] into BriK^ — > 7r°7r£,) = 2xBr(Ki — > 
7r°7r°) x Br(7T° — > 77) x Br(n° — > e + e~7), the single event sensitivity was ob- 
tained to be (3.03 ± 0.43) x 10~ 7 , where the error represents the statistical 
uncertainty. The upper limit on the branching ratio was determined to be 

Br(K L -> 7r°7e + e") < 7.1 x 10~ 7 (90% C.L.), 



in which the statistical error of the normalization mode was taken into ac- 
count [14]. 

In summary, we performed an experimental search for the Kl — > 7r°7e + e~ 
mode. We observed no events and set a 90% confidence level upper limit of 
Bt(Kl — > 7r°7e + e~) < 7.1 x 10~ 7 for its branching ratio. This is the first 
published experimental result on this decay mode. 
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